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Abstract 

We present a model for lepton masses and mixing angles using the deconstruction 
setup, based on a non-Abelian flavor symmetry which is a split extension of the 
Klein group Z2 x Z2. The symmetries enforce an approximately maximal atmo- 
spheric mixing angle 623, a nearly vanishing reactor mixing angle ^13, and the hier- 
archy of charged lepton masses. The charged lepton mass spectrum emerges from the 
Froggatt-Nielsen mechanism interpreted in terms of deconstructed extra dimensions 
compactified on 5^. A normal neutrino mass hierarchy arises from the coupling to 
right-handed neutrinos propagating in latticized /Z2 orbifold extra dimensions. 
Here, the solar mass squared difference Auiq is suppressed by the size of the dy- 
namically generated bulk manifold. Without tuning of parameters, the model yields 
the solar mixing angle 612 = arctan l/V'2'- Thus, the obtained neutrino masses and 
mixing angles are all in agreement with the Mikheyev-Smirnov-Wolfenstein large 
mixing angle solution of the solar neutrino problem. 
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1 Introduction 



The standard model of elementary particle physics (SM), minimally extended 
by massive Majorana neutrinos, contains 22 fermion mass and mixing param- 
eters (6 quark masses, 6 lepton masses, 3 CKM mixing angles [1,2], 3 MNS 
mixing angles [3], 2 Dirac CV violation phases, and 2 Majorana phases). In 
four- dimensional (4D) gauge theories the fermion masses and mixing angles 
result from Yukawa couplings between fermions and scalars. However, even 
in 4D grand unified theories (GUTs) the Yukawa couphngs are described by 
many free parameters and one therefore often assumes that the structure of 
the Yukawa coupling matrices is dictated by a sequentially broken flavor sym- 
metry. Indeed, since the 4D Yukawa couplings become calculable in higher 
dimensional gauge theories [4], the effective fermion mass matrices emerging 
after dimensional reduction should actually be highly predictable from sym- 
metries or the topology of internal space [5] . 

At first sight, the hierarchical structure of CKM angles and charged fermion 
masses suggests an underlying non-Abelian flavor symmetry group which es- 
sentially acts only on the flrst and second generations [6]. However, in these 
theories it is difficult^ to obtain exact predictions compatible with recent 
atmospheric [8,9] and solar [10,11] neutrino experimental results which pre- 
fer the Mikheyev-Smirnov-Wolfenstein (MSW) [12] large mixing angle (LMA) 
solution of the solar neutrino problem [13]. In fact, the KamLAND [14] reac- 
tor neutrino experiment has recently conflrmed the MSW LMA solution at a 
signiflcant level [15,16,17,18]. In the "standard" parameterization, the MSW 
LMA solution tells us that the leptons obey a bilarge mixing in which the solar 
mixing angle 612 is large, but not close to maximal, the atmospheric mixing 
angle 623 is close to maximal, and the reactor mixing angle 6*13 is small. 

By assuming only Abelian U{1) [19] or Z„ [20] symmetries one flnds that the 
atmospheric mixing angle may be large but cannot be enforced to be close to 
maximal. Thus, a natural close to maximal i/^-i/^-mixing can be interpreted 
as a strong hint for some non-Abelian flavor symmetry acting on the 2nd and 
3rd generations [21,22,23]. Neutrino mass models which give large or maximal 
solar and atmospheric mixing angles by assigning the 2nd and 3rd generations 
quantum numbers of the symmetric groups S2 [24] or S3 [25] have, in general, 
difficulties to address the hierarchy of the charged lepton masses. This problem 
can be resolved in a supersymmetric model for degenerate neutrino masses 
by imposing the group A4, the symmetry group of the tetrahedron [26]. In 
addition, this model predicts exactly 623 = f from the A4 symmetry and 
gives with some tuning of parameters the solar mixing angle 612 of the MSW 
LMA solution. However, in unified field theories a normal hierarchical neutrino 



^ For a possible counter-example, see the recent model in Ref. [7]. 
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mass spectrum seems more plausible than an inverted or degenerate one [27]. 
A comparably simple way of accommodating the MSW LMA solution for 
normal hierarchical neutrino masses is provided in scenarios of single right- 
handed neutrino dominance [28]. 

Although unification in more than four dimensions can serve as a motivation 
for flavor symmetries, higher-dimensional gauge theories have dimensionful 
gauge couplings and are usually non-renormalizable. They require a trunca- 
tion of Kaluza-Klein (KK) modes near some cut-off scale Mf at which the 
perturbative regime of the higher-dimensional gauge theory breaks down. Re- 
cently, however, a new class of 4D gauge-invariant field theories for decon- 
structed or latticized extra dimensions has been proposed, which generate the 
physics of extra dimensions in their infrared limit [29,30]. These theories are 
renormalizable and can thus be viewed as viable UV completions of some fun- 
damental non-perturbative field theory, like string theory. In deconstruction, 
a number of replicated gauge groups is connected by fermionic or bosonic link 
variables. Since replicated gauge groups also frequently appear in models using 
the Froggatt-Nielsen mechanism [31,32], it is straightforward to account for 
the generation of fermion mass matrices by deconstructed extra dimensional 
gauge symmetries. 

In two previous works [33,34] we have presented versions of a model for bi- 
large leptonic mixing based on a vacuum alignment mechanism for a strictly 
hierarchical charged lepton mass spectrum and an inverted neutrino mass hi- 
erarchy. As a result of the inverse hierarchical neutrino mass spectrum, we 
obtained the approximate relation ^12 = f — ^13 and hence the lower bound 
37° < 9i2 on the solar mixing angle. In this paper, we now apply this vacuum 
alignment mechanism to a model for lepton masses and mixing angles which 
yields more comfortably the MSW LMA solution with normal hierarchical 
neutrino masses. Specifically, the lepton mass hierarchies are generated by de- 
constructed extra-dimensional U{1) gauge symmetries where the lattice-link 
variables are themselves subject to a discrete non-Abelian flavor symmetry. 

The paper is organized as follows: in Section 2, we introduce our example field 
theory by briefiy reviewing the periodic and the aliphatic model for decon- 
structed extra dimensions. In Section 3, we first assign the horizontal U{1) 
charges, then we discuss the non-Abelian discrete flavor symmetry and exam- 
ine its normal structure. Next, in Section 4, we construct the potential for the 
scalar link- variables of the latticized extra dimensions from the representations 
of the dihedral group ^4. The vacuum structure which results from minimizing 
the scalar potential is determined and analyzed in Section 5. In Section 6, we 
describe the generation of the charged lepton masses via the vacuum alignment 
mechanism. Then, in Section 7, we study the neutrino mass matrix, determine 
the neutrino mass and mixing parameters and match the types of latticiza- 
tions of the orbifold extra dimensions onto the presently allowed ranges for 
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Fig. 1. Moose diagram for the gauge group = Ilf^iGj. 




Fig. 2. Moose diagram for the gauge group = Ilf^iG^. 

AttIq impUed by recent KamLAND results. Finally, we present in Section 8 
our summary and conclusions. In addition, Appendix A gives a brief review 
of the dihedral group ^4 and in Appendix B the minimization of the scalar 
potential is explicitly carried out. 



2 Deconstruction setups 

2.1 The periodic model 

Consider a 4D gauge-invariant field theory for deconstructed or latticized ex- 
tra dimensions [29,30]. Let the field theory for j = 1,2, .. . be described by 
products = nZ^Gl of U{1) gauge groups Gj , where Nj is the total number 
of "sites" corresponding to G^ . In the "periodic model" [29], one associates 
with each pair of groups x Gj^_^_l a link-Higgs field Qj with charge (-1-1, — 1) 
under G^ x G^+i, where i is periodically identified with i + Nj. For our example 
field theory we will assume one periodic model for each of the gauge groups 
G^ and G^ where the number of sites is iVi = 2 and N2 = 4. This field theory 
is conveniently represented by the "moose" [35] or "quiver" [36] diagrams in 
Figs. 1 and 2. Restricting for the present to a single product gauge group G^ 
with Nj — N sites, we can drop the index j and write the Lagrangian in the 
periodic model as 

1 Af N 

=^ = -T E Fif^'^F'"' + E i^^Qi)^ D^Qi^ (1) 
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di,Ai^ is the 4D field strength and the covariant derivative 



where F^^,, 
is defined by 

Df,Qi = (d^ - igiAi^ + igi+iA(^i+i)^)Qi. (2) 
where gi is the dimensionless couphng constant of the gauge group G at the 
ith site. For simplicity, it is assumed that all gauge couplings are equal, 



91 = 92 



9n = 9- 



(3) 



When the Qi aquire after spontaneous symmetry breaking (SSB) universal 
vacuum expectation values (VEVs) and each field becomes a non-linear a- 
model field 

Qi ^ V exp{iai/v), (4) 
the full gauge symmetry is higgsed to the diagonal U{1) and it is seen that 
the Lagrangian in Eq. (1) describes a 4+1 dimensional U{1) gauge theory 
on a flat background, where only the fifth dimension has been latticized. The 
lattice spacing a and circumference R of the fifth dimension are a = 1/ (gv) and 
R = N/ (gv), whereas the 5D gauge coupling g^ is given by = 1/ {ag'^) [29, 
37]. (In the model with free boundary conditions [30] generic gauge couplings 
and non-universal VEVs introduce an overall non-trivial warp factor [38,39], 
resulting, e.g., in a Randall- Sundrum model [40].) Implicitly, we identify the 
link-Higgs fields with the continuum Wilson lines 



f.(i+l)a 



exp 



(5) 



where {x^,x^^) are the bulk coordinates and is the fifth component of the 
bulk gauge field associated with the gauge group G. The gauge boson mass 
matrix takes the form 



N 



(6) 



i=l 



and is of the type of a nearest neighbor coupled oscillator Hamiltonian. In the 
case of iV = 4, for example, the mass matrix for the gauge bosons reads 



2 2 
9 V 



^2-10 -1^ 
-12-10 

0-12-1 
y-1 0-12 



(7) 



In general [29], the gauge boson mass matrix yields a mass spectrum labeled 
by an integer k satisfying —N/2 < k < N/2, 



Ml = Ag'^v'^sm^ 



sm 



(8) 
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Fig. 3. Moose diagram for bulk vector fields and fermions in the aliphatic model, 
where ps = 27rk/R is the discrete 5D momentum. For small |A;| the masses are 

277-1 i; I 

Mk ^ \Pk\ = 1^1 « N/2, (9) 

which is exactly the KK spectrum ^ of a 5D gauge theory compactified on 
with circumference R = N/ (gv). The zero mode corresponds to the unbroken, 
diagonal U{1). Hence, at energies E <^ gv / N we observe an ordinary 4D gauge 
theory, in the range gv/N < E < gv the physics is that of an extra dimen- 
sion, and for E ^ gv an unbroken U{1)'^ gauge theory in four dimensions is 
recovered. 



2.2 The aliphatic model for fermions 



The "aliphatic model" [30] for some product gauge group {j = 3,4, . . .), 
as defined in Sec. 2.1, is obtained from the periodic model by setting Q-'^^. = 
which yields a linear system with free boundary conditions. For this type of 
latticization we will consider Nj SM singlet Dirac fermions ^{(n = 1, . . . , A^^), 
each of which carries an associated charge —1. The SM singlet Higgs fields 
Qn = 1, . . . Nj) which are assigned the x Gl^^i charges (+1, —1) specify 
the allowed couplings between the fermions. Restricting here to a single prod- 
uct gauge group G^ we may in this section drop the index j for the rest of the 
discussion. The moose diagram for the ahphatic model is shown in Fig. 3. We 
will denote by and "^nR the left- and right-handed chiral components of 
respectively. To engineer chiral fermion zero modes from compactification 
of the 5th dimension, one can impose discretized versions of the Neumann and 
Dirichlet boundary conditions — '^{n-i)l — and = 'iNR = ^ which 
explicitly break the Lorentz group in five dimensions. Using the transverse 
lattice technique [41,42] the relevant mass and mixing terms of "ifnL and "ifmR 
then read [30, 38] 



Af-l 
n=2 




^(n+l)i? - *ni?) - ^nR (^nL - ^*(n-l)L 



+Mf^lL'^2R + ii-C., (10) 



Note here the doubling of KK modes. 
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where v = (Qn) for n = 1, . . . , — 1, i.e., we have universal VEVs. In fact, 
Eq. (10) is the Wilson-Dirac action [43] for a transverse extra dimension which 
reproduces the 5D continuum theory in the hmit of vanishing lattice spacing. 
Clearly, this tacitly presupposes a specific functional measure for the link 
variables in question which may, however, no longer be a necessary choice 
when the lattice spacing is finite [42] . After SSB the mixed mass terms of the 
chiral fermions are 



N-l 



^mass = M/ [^ni(*(n+l)i? " *nil) " *nil(*nL - in-l)L] 



n=2 



= (^iL, . . . , ^(iv-i)L)^A^i(*2ii, ■ ■ ■ , ^{n-i)r) + h-c, 
where the (A^ — l)x(A^ — 2) fermion mass matrix M.i is on the form 



/ 



1 
-1 1 

00 







-1 1 

-1 



fill 



(12) 



Diagonahzing the {N — 2) x {N — 2) matrix MiMi gives for the mass eigen- 
values of the right-handed fermions 



Mr,R = 2Mfsmi ^^_ \, n=l,2,...,7V-2, (13) 
where the associated mass eigenstates '^kR are related to by [30]: 



^nR = Y E [^J^j "^kR. (14) 



k=l 



The diagonahzation of the (A^ — 1) x (A^ — 1) matrix MiM-l yields masses for 
the left-handed fermions, which are identical with the gauge boson masses, 

2^N_l) ]^ n = 0,l,...,Ar-2, (15) 

and in terms of the associated the mass eigenstates '^kL we have [30] : 



/AT - 1 ^2 /2n + l A;7r \~ 



fc=0 
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Setting Mf — and taking the limit n <^ N the hnear KK spectrum of a 5D 
fermion in an orbifold extra dimension S^/Z2 is reproduced. Up to an overall 
scale factor of 2, the periodic and the aliphatic model generate identical KK 
mass spectra for bulk vector fields, with the number of KK modes doubled in 
the periodic case [44] . 



3 Horizontal charges 

In 5D continuum theories hierarchical Yukawa coupling matrix textures have 
been successfully generated [45] . Hierarchical Yukawa matrices have also been 
obtained in deconstructed warped geometries [46] . In our deconstruction setup 
of Sec. 2 let us now consider an extension of the SM, where the lepton masses 
arise from higher-dimensional operators [47], partly via the Proggatt-Nielsen 
mechanism. 



3.1 U{1) charges 

We denote the left-handed SM lepton doublets as La = {vaL Gcxl)^ and the 
right-handed charged leptons as Ea, where a = e,/x, r is the flavor index. 
For simplicity, we will assume the electroweak scalar sector to consist only 
of the SM Higgs doublet H. In order to account for the neutrino masses, 
we will additionally assume three SM singlet scalar fields $,o,Ci, and {2, as 
well as three heavy SM singlet Dirac neutrinos Fe,Ff^, and Fr. Since these 
Dirac neutrinos are supposed to have masses of the order of the GUT scale 
~ 10^^ GeV, they will account for the smallness of the neutrino masses via 
the seesaw mechanism [48,49]. Furthermore, the standard Froggatt-Nielsen 
mechanism is implied in the charged lepton sector by the presence of heavy 
fundamental charged fermion messengers. The electron, muon, and tau masses 
will be denoted by me,m^, and m,-, respectively. 

In models of inverted neutrino mass hierarchy, a small reactor mixing angle ^13 
can be understood in terms of a softly broken L^ — L^^ — Lr lepton number [50]. 
Analogously, we assume that our example field theory contains a product 

X Gr of two extra U{1) gauge symmetries Gl and Gr which distinguish 
the 1st generation from the 2nd and 3rd generations, but act also on Ci) ^2, 
and the scalar link fields of the product gauge groups G^ — 11^^^ and 

= Tlj^^Gf. The x charge assignment is shown in table 1. Note that 
the Gl x Gr symmetry is anomalous. However, anomalous U{1) charges are 
expected to arise in string theory and must be canceled by the Green-Schwarz 
mechanism [51]. 
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GlxGr 



Ee, Fe, Co 

Lh,Lt,Q\,QI, ^1,^2 

F F 

^ fit 



QlQl 



(+2,-2) 
(-2, +2) 
(+1,0) 
(-1,0) 

(o,+i) 
(o.-i) 



Tabic 1 

Assignment of the horizontal U{1)'^ charges to the different fields. The fields not 
shown here transform trivially under Gl x Gr. 




(H) {QD {QD 




/ \ \ 



{QD 



{QD 



Fig. 4. Lowest-dimensional Yukawa interactions generating for a,/3 
charged lepton mass terms in the /n-r-subsector. 



{QD 



II, T the 



Let us now suppose that the charges of the product gauge groups Gl x Gr, G^, 
and G^ are approximately conserved in the charged lepton sector, implying the 
generation of hierarchical charged lepton mass terms via the Froggatt-Nielsen 
mechanism. Furthermore, we assume that the fields ^Oi 'Cii ^"^^ ^2 carry nonzero 
G^ and G^ quantum numbers , respectively, while the fundamental charged 
Froggatt-Nielsen messengers transform only trivially under G^ and G^. As a 
consequence, the fields Co,Ci, S can be discarded in the discussion of 
the generation of the charged lepton masses (sec also Sec. 7.1). Then, we 
conclude from table 1 that gauge-invariancc under the product gauge group 
Gl X Gr X G^ X G'^ allows in the /i-r-subscctor only two general types of 
non-renormalizable leading-order charged lepton mass terms: One operator 
of the dimension 4: + Ni — 6 and one of the dimension 4 + A^2 = 8. The 
corresponding Froggatt-Nielsen-type diagrams are shown in Fig. 4. Note here, 
that the effective Yukawa couplings of the dimension-six and dimension-eight 
operators arise from the link fields of the deconstructed extra-dimensional 
gauge groups G^ and G^, respectively. 



The exact G^ and G^ charge assignment to these fields is discussed in Sec. 7.1. 
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3.2 Discrete charges 



It has been pointed out, that a naturally maximal i/^-i/^-mixing is a strong 
hint for an underlying non-Abelian flavor symmetry acting on the 2nd and 
3rd generations [21,22,23]. Seemingly, this flavor symmetry is unlikely to 
be a continuous global symmetry because of the absence of such symmetries 
in string theory [52]. From the point of view of string theory, however, it 
is interesting to analyze conventional grand unification in connection with 
discrete symmetries since these might provide a solution to the doublet-triplet 
splitting problem [53,54,55]. Additionally, discrete symmetries can ensure an 
approximately fiat potential for the GUT breaking fields in order to avoid 
fast proton decay [56]. Moreover, it has recently been demonstrated that the 
appropriate discrete symmetries of the low-energy theory can be naturally 
obtained in deconstructed models with very coarse latticization, when the 
lattice link variables are themselves subject to a discrete symmetry [55]. It 
is interesting to note, that although a classical lattice gauge theory based 
on a system of link fields which are elements of some discrete group has no 
continuum limit, this does not necessarily carry over to the quantum theory^ . 

Among the discrete symmetries that have been proposed in the context of 
the MSW LMA solution as a possible origin of an approximately maximal 
atmospheric mixing angle 6*23 are the symmetric groups S2 and S3 acting on 
the 2nd and 3rd generations of leptons [25,24]. While this can indeed give an 
approximately maximal f^-i^r-mixing, the hierarchical charged lepton mass 
spectrum is typically not produced in these types of models, since they rather 
yield masses of the muon and tau that arc of the same order of magnitude. 
However, when we add to the regular representation of S2 in the /x-r-fiavor 
basis the generator diag(— 1, 1), one obtains the vector representation of the 
dihedral group ^4, which is non-Abehan (see App. A). Since the generator 
diag(— 1, 1) distinguishes between and Lr it may serve as a possible source 
for the charged lepton mass hierarchy, which breaks the permutation symmetry 
S2 C ^4 characteristic for the i/ju-^',--sector. Clearly, if the charged lepton 
masses arise from the Froggatt-Nielsen mechanism, we will have to expect 
that the underlying symmetry is actually equivalent to a group extension^ 
of some subgroup of ^4, presumably a suitable subgroup of some replicated 
product of ^4-factors. 

Motivated by these general observations, we will take here the view, that 
an approximately maximal atmospheric mixing angle |^23 — f | <S 1 is due 
to a non-Abelian discrete symmetry W which is a group extension based on 



^ For a pedagogical introduction to lattice gauge theory see, e.g., Ref. [57]. 

^ An extension of a group by a group H is an embedding of N into some group 

E such that N is normal in E and H ~ E/N. 
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the generators of the dihedral group ^4 under which the leptons of the 2nd 
and 3rd generation, the scalars ^i,^, and the scalar link variables of the 
deconstructed extra-dimensional gauge symmetries = U^^iG] and ~ 
Uf^^Gj transform non-trivially. Specifically, we suppose that in the defining 
representation V the group which is a subgroup of the four-fold (external) 
product group ^4 x ^4 x ^4 x ^4, can be written as a sequence 



V{g) = diag 



D^{g), D^ig), D^{g), D,{g) 



9& 



(17) 



where each element g E ^ is associated with four^ (in general different) 
2x2 charge operators Di{g) (i — 1,2,3,4) which have values in the vector 
representation of ^4. Clearly, for given i = 1,2,3,4 the set of all operators 
{Di{g) : g G } forms a subreprescntation of which we will call Di. Now, 
using the notation of App. A we suppose that ^ is represented by four gener- 
ators with following ^4-charge structure 



V{ai) = diag 



V{a2) =diag 



V{a3) =diag 



V^a^) = diag 



D{C,), D{E), D{C,,), D{C,, 



D{Cv), D{C,,), D{Cb), D{E) 



D{E), D{E), D{C,), D{C,) 



D{C,), D{C,), D{E), D{E) 



(18a) 
(18b) 
(18c) 

(18d) 



where 01,02,03,04 G ^ are the corresponding abstract generators. Note that 
these operators are characterized by a one-to-one-correspondence 'P(oi) ^ 
T>{a2) and 1^(03) T>{a^ under the permutation of the upper-left and the 
lower-right 4 x 4-matrices. As will be shown in Sec. 3.3, by factoring the 
subreprescntation Di for any i = 1,2,3,4 with respect to its kernel A/i we 
obtain a representation of the factor group ^/J\fi that is isomorphic with ^4. 
This implies, of course, that all four subrepresentations Di, D2, -D3, and D4 are 
two-dimensional irreducible representations (irreps) of ^. From App. A it is 
seen, that in an appropriate basis the generators of Eqs. (18) can be explicitly 
written as 



^ Since each element g ^ is associated with four (in general) different operators 
Di{g), D2{g), D3{g), and D^i^g), one may view in a discretized picture P as a 4- 
valued representation of ^4 with ^ as the corresponding (universal) covering group. 
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'D(ai)= dia; 



V{a2) =diag 



D(a3) =diag 



V( 




1 


(::) 








(io\ 




,0 ij 




(19a) 



(19b) 



(19c) 



(19d) 



With respect to V we will combine the left- and right-handed SM leptons 
as well as the Dirac neutrinos of the 2nd and 3rd generations into the ^- 
doublets 2^ = (L^, L^)^, 2^ = {E^, E^Y , and 2f = {Fi^, F^Y, respectively. 
The fermionic doublets 2^ and 2^;^ as well as the scalars = (^la, Cib)'^ and 
^2 = (^^20! '^26)"^ arc all put into the doublet representation Di. Next, the gen- 
eralized Wigner-Eckart theorem tells us that the effective Yukawa interaction 
matrix spanned by 2^ and 2^ in the /^-r-subsector of the charged leptons is 
identical with a linear combination of sets of irreducible Yukawa tensor opera- 
tors. If the irreducible Yukawa tensor operators take their values in the vector 
representation of ^4, it is clear from App. A, that the hierarchy <C m-r is 
only possible if 2^; transforms according to an irrep of ^ which is inequiva- 
lent ^ with Di. We will therefore put 2e into the irrep D2 which is inequivalent 
with Di and note that the large subgroup of ^ generated by 1^(02), ^^(03), and 
X>(a4) acts diagonally on 2^ and 2^. 

In Sec. 3.1 we have seen that GlX Gr xG^ xG'^ gauge- invariance requires the 
possible lowest-dimensional effective Yukawa interactions between 2^ and 2^ 
to be identical with the dimension-six and dimension-eight Proggatt-Nielsen- 
type operators displayed in Fig. 4. In theory space, we can therefore identify 
the sets of irreducible Yukawa tensor operators spanned by 2^ and 2e with the 
Wilson loops around the plaquettes associated with the deconstructed extra- 
dimensional gauge-symmetries G^ and G^: The dimension-six and dimension- 
eight operators correspond to the Wilson loops around the moose diagrams 
of G^ (Fig. 5) and G^ (Fig. 6), respectively. Although the irreps Di and D2 
already determine the overall transformation properties of the Wilson loops 
under there is still some ambiguity in the individual ^-charge assignment 



^ Similarity transformations allow only mappings within one class, which would 
yield = rrij- after SSB. Note also that different transformation properties of left- 
handed and right-handed fermions under horizontal symmetries are used in models 
of "neutrino democracy" [64]. 



12 




Fig. 5. Contour in theory space which generates the dimension-six mass operator 
in Fig. 4 by connecting 2^ and 2e via the hnk fields of (left panel). Formal 

contraction of the open fermion lines and Gl x Gr into G\ (dashed boxes) exhibits 
the full G^ gauge-invariance of the corresponding effective Yukawa couplings (right 
panel). 




Fig. 6. Contour in theory space which generates the dimension-eight mass operator 
in Fig. 4 by connecting 2^ and 2e via the link fields of (left panel). Formal 
contraction of the open fermion lines and Gl^ Gr into G\ (dashed boxes) exhibits 
the full G^ gauge-invariance of the corresponding effective Yukawa couplings (right 
panel). 

to the involved link-fields Q\^Q\., and Q\^Q\^Q\^Q\- Here, it is appealing 
to assume that all the scalar link fields of and transform according 
to some doublet subrepresentation of which has the property that its lift 
is isomorphic with ^4. Furthermore, a non-trivial Yukawa matrix structure 
requires that the products 11^^^ and 11^^^ are non-singlet representations 
of ^. Since each of the Wilson loops involves an even number of (two or 
four) link fields, we find from the multiplication rules in App. A that this can 
only be the case if each of the sets {Q\iQ:\} and {Q\-,Q2-,Q\-iQ^a} contains 
at least two inequivalent irreps. This can be simply realized by putting, e.g., 
the 1st (i = 1) link fields Q\ and Q\ into the doublet representation while 
the remaining link fields Q\,Q\,Q\, and Q\ all transform as doublets under 
1^4. Note again, that both and -D4 are characterized by a large common 
subgroup generated by ©(ai), ©(03), and 1^(04), which acts diagonally on all 
of the link fields. In component-form, the scalar ^-doublets will be written 
as Q\ = {qi, qlbf, where i = 1,2, and = {q^, qlf, where i = 1,2,3,4. 
To complete the ^-charge assignment, we suppose that the first generation 
fields Lg, Eg, and F^, as well as the scalar ,^0 transform only trivially under $f , 
i.e., they are all put into the identity representation / of ^. For the fermionic 
^-singlets Lg, Ee, and Fg we will also choose the notation Ig = -Lg, 1b = -E'e, 
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Di 


D2 


D3 


L>4 


le, 1e, 1f,?0 


2^, 2f,^i,^2 


2e 


Q\, Qi 


Q2! Qi 



Table 2 

Assignment of the fermionic and scalar fields to the different irreps and /. Note 
that the scalar link fields transform according to their position in the "index space" 
of gauge groups. 

and Ip = Fe- The assignment of the fields to the irreps {i — 1, 2, 3, 4) and 
/ is summarized in table 2. 

In this subsection, we have presented the non-Abelian group in terms of its 
generators which are motivated by phenomenology. In the following subsection 
we will examine in some more detail the structure of ?f through a descending 
series of normal subgroups, i.e., the normal structure. 

3.3 Normal structure 

In Sec. 3.2 the discrete group ^ has been constructed from the generators of 
the dihedral group ^4. A standard way of gaining further information about 
^ is to analyze a series of subgroups of $f , where each term is either normal 
in ^ or at least normal in the previous term. In general, if a subgroup of 
some group G is normal in G we will write N < G. 

Wc will denote by /Ci the collection of all elements g & ^ which obey Di{g) = 
I2, where I2 is the 2x2 unit matrix. Accordingly, we will define /C2 as the 
set of all elements g & ^ which obey Di(g) = D2{g) = I2, and /C3 as the set 
of all elements g & which obey Di{g) — D2{g) — D3{g) — I2. We therefore 
have the sequence 

/C3 C /C2 C /Ci C ^, (20) 

where each subset is a group, actually an invariant subgroup of the embedding 
groups, for if a G /Cj (i = 1, 2, 3) is homomorphically mapped on the identity 
of the operator group Di (B . . . Q) D^, then so are all elements in its class. 
We therefore have /Cj+i < JCi and /Cj < ^ for every i, i.e., the subgroup 
series in Eq. (20) is in fact a normal series. Prom Eqs. (19) we find that for 
any k E IC2 the decomposition of T>{k) into a product of the operators T>{ai) 
{i = 1,2,3,4) necessarily involves each of the factors V{ai),V{a2), and ©(04) 
an even number of times. In turn, this implies that the operators D^ik) and 
Di{k) are on diagonal form and can take their values only in the classes C|, 
and 6*2(2) of the dihedral group ^4 (see App. A). Since the number of elements 
in the classes Ci, and 6*2(2) is four, we conclude that the order of /C2, which 
we win denote by I/C2I, obeys I/C2I < 4. Indeed, besides the unity, we find the 
three distinct elements and {ttia^Y which are all contained in K,2. 

Hence, I/C2I =4 and we have 
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/C2 ~ {diag [D{E), D{E), D{E), D{E)] , 



di^g[D{E),D{E),D{Ca),D{Ca)], 



die.g[D{E),D{E),D{C,),D{C,)], 



diag[D{E),D{E),D{Cl),D{Cl)]}. 



(21) 



Since /C3 C IC2 the group /C3 is trivial and contains only the identity of V, 
as inspection of Eq. (21) shows. From Eq. (21) we find that IC2 describes a 
2- fold axis with a system of two 2- fold axes at right angle to it. Hence, JC2 is 
isomorphic with the Klein group Z2X Z2 — /C2. In fact, the Klein group is one 
of the dihedral groups Z2 x Z2 ~ ^2 (see App. A) . With /C2 exphcitly given 
in Eq. (21) we can now easily construct a representation of the embedding 
group /Ci ^ /C2- First, note that for any k E Ki the resolution of the operator 
T){k) into a product of the generators in Eq. (19) must involve 'D{a2) an even 
number of times. Hence, D2{k) is necessarily on diagonal form implying that 
the index of the subgroup /C2 under the group /Ci is at most four. In fact, the 
group /Ci can be decomposed into (right) cosets ^ in terms of 



where we can choose the different coset representatives to be hi — 0104, 62 — 
02^102, and 63 = 62^1- In other words, four elements of /Ci are mapped on 
each element of the operator group corresponding to the representation of 
/Ci subduced by D2 (subduced representation) D2 [ ICi. Since the subgroup 
K,2 ^ /Ci is the kernel of -D2 i ^1, the operators D2{bi) {i — 1, 2, 3) from the 
(right) transversal allow us to identify the lift of D2 i ]Ci sts the Klein group 
Z2 X Z2. As for the coset representatives don't commute with the elements 
of /C2 it follows that /Ci is actually an (external) semi-direct product or split 
extension of JC2 by Z2 x Z2. Hence, we can write 



where ip is the mapping of IC2 into the automorphism group of IC2 which is 
determined by the particular choice of the generators in Eq. (18). The mapping 
(f can be identified with the conjugation homomorphism of the corresponding 
internal semi-direct product describing the interaction of the two involved 
Z2 X Z2 subgroups inside ^. Now, the order of /Ci is simply given by the 
product of the orders of the semi-direct factors, i.e., |/Ci| = 16. Since we have 
^//Ci ~ ^4 the decomposition of into right cosets with respect to JCi reads 



^ Since the subgroup /C2 is invariant, left and right cosets are identical. 
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/Ci~/C2 + E^2&, 



(22) 



/Ci ~ {Z2 X Z2) x^ {Z2 X Z2), 



(23) 
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^~/Ci + X:/CiQ, (24) 

where we can choose for the coset representatives Ci = ai, C2 = 020102, C3 = 
(0102)^, C4 = 02, C5 = 0102, Cq = 0201, and C7 = 010201. Here, the irrep Di 
maps 16 elements of onto each element of ^4 and hence |^| = 128. Again, 
the seven coset representatives don't commute with the elements of /Ci which 
shows that $f is a split extension of /Ci by ^4, i.e., 

^ ~ (Z2 X Z2) (Z2 X Z2) x^ ^4, (25) 

where ■?/' is identified with the conjugation homomorphism describing the in- 
teraction of ^4 with the involved Z2 x Z2 subgroups inside By construction, 
the operator group associated with the irrep Di has ICi as its kernel. Hence, 
by factoring Di with respect to the subgroup /Ci we see that the lift of Di 
is isomorphic with ^4. Furthermore, by replacing X'(oi) — * I'(oi)X>(a4), we 
confirm that the irreps Di and D2 change their roles in the above consider- 
ations. Taking, in addition, the exchange symmetry between 1^(01) 1^(02) 
and 1^(03) ^ 1^(04) under permutation of the 4x4 submatrices on the diag- 
onal into account (see Sec. 3.2), one generally concludes that the hft of every 
subrepresentation Dj (i = 1, 2, 3, 4) is isomorphic with ^4. 

In this subsection we have analyzed the structure of the group ^ and its 
relation to the dihedral groups ^2 and ^4. Specifically, we have seen that the 
lift of any irrep (i = 1, 2, 3, 4) is isomorphic with ^4. This will allow us in 
the next section to apply the decomposition and multiplication rules of ^4 in 
order to find the relevant ^-singlets in product representations. 



4 Construction of the scalar potential 

We will denote by $ and Q two arbitrary scalar $f -doublets which are listed 
in table 2 and write them in component-form as 

$ = (0„, (f),f , n = {Ua, CUhf , (26) 

where $, f2 e {^1,^2, QliQhQiiQhQliQl}- definiteness, we will assume 
that $ is put into the irrep and O is put into the irrep Dj, where — 
1, 3, 4. Gauge-invariance under the C/(l) product groups ]Tj^^G^ and Gl x Gr 
tells us that any renormalizable term in the scalar potential which involves, 
e.g., the field $ must actually contain the tensor product $ (S> The lowest- 
dimensional ^^-invariant operator-products of $ in the multi scalar potential 
are therefore 

Vo(^) = $ $t|^^ ^ Ao (|(/>a|' + l^fel') , (27) 
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where $ transforms under any of the irreps Di {i — 1,3,4) and Aq is a real- 
valued number. Note that three-fold products of $ and in the potential are 
forbidden by both the U{1) charge assignment as well as by ^^-invariance. Ac- 
cordingly, any invariant term in the scalar potential which mixes the doublets 
$ and Q must be contained in the product representation $(8)$^(8)Q(8)Q^ Since 
the liftings of the irreps Di and Dj are isomorphic with ^4 (see Sec. 3.3), the 
^-invariant mixed terms of $ and VL are found by considering all combinations 
of the one- dimensional irreps of ^4 in the product 

$ ® $t ^ Q ® ^ [li + 1^ + l^c + I'd] ® [1a + 1b + Ic + Id] > (28) 

where each of the sets IajIr' Idj ^"^^ IajIb' Ic'Id respectively denotes 
the ^4-singlet representations associated with Di and Dj (see App. A). In 
component form, the singlet representations are explicitly given in Eq. (A. 3). 
In order to extract from Eq. (28) the relevant dimension-four terms we will 
first suppose that $ and Q, are put into the same irrep D; = Dj. Then, the 
decomposition of the product representations in Eq. (A. 2) in conjunction with 
the multiplication rules in table A. 2 yield that the invariant mixed terms of 
$ and fl are in this case 



Vi($,l]) = ($®$t)^(^^^t)|^^ 

= A (|0ar + |<^6|') (|c^„|' + |c^6r) 
+Bi (^(l)a(l)l - (t>l(t>b) [^a^^l - i^li^b) 

+C4 - (Kr-|a;,r) 

+ ((>%) (ujaujl + ulub) , (29a) 

where we have labeled the real- valued coefficients Ai, Bi,Ci, and Di of the in- 
variants according to the sequence of products (1a)^i (1b)^i (lc)^i and (1d)^- 
Let us now turn to the case, when $ and ^2 belong to different irreps D^ ^ Dj. 
First, suppose that $ transforms under D3 and Q transforms under D4. Ap- 
phcation of the generators I>(ai), I>(a2), and X>(a3) yields that in Eq. (28) for 
i = 3 and j = 4 only the combinations 1^ (E) 1^ and 1^ ® Ic are ^-invariants. 
Hence, the most general invariant mixed term of $ and Q reads 



V,(^, n)=A, + + 

+C2 (|0„r - \<Pb\') (kal' - , (29b) 

where $ G {Q\, Qj}, ^ G {Ql, Ql, Q|, Ql} and the coefficients A2, C2 are real- 
valued numbers. Now, suppose that $ transforms under Di and Q transforms 
under one of the irreps D3 or D4. In this case, subsequent application of the 
operators T>{ai) {i — 1, 2, 3, 4) to Eq. (28) readily yields that the most general 
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mixed term of the fields $ and Q is given by 1a ® 1a hence 

Vs{^, n) = As {\<Pa\' + l^^r) (kar + , (29c) 

where $ e {^1)^2}, ^ € {Qi, Qi) Qi) Qi; Ql} ^-^d the coefiicient A3 is 
a real-valued number. Putting everything together, the most general multi 
scalar potential V of the SM singlet scalar fields decomposes into the terms 
Vo{^)yi{^,n), V2{^,n), and Vsi^,^) of Eqs. (29), as follows 



V=J2[Vo(^,) + Vo(Q]) + Vo(Ql)]+ E Vi(^h,U 

i,j,k 11,12=1,2 

+ Yl V^{X,Y)+ [Vl{Ql,Ql) + V^{Ql,Q 



2 

4 



X,y=Q},Q2 jij2=2 



+vi(Q2, Q2) + E E [^2(x, Q^) + 1-2(^, g|i)" 

X=Qi,Qf ji=2 

+ E[^3(6,g]) + ^3(6,gD], (30) 

where i = 1,2; j' = 1,2, and k — 1,2,3,4. In Eq. (30) we have omitted 
and the SU (2) Higgs doublet H. Actually, in any renormahzable terms of the 
full multi-scalar potential which mix the ^-singlets with the ^-doublets, the 
^-singlet fields .^o H are only allowed to appear in terms of their absolute 
squares |^oP and \H\'^. This is an immediate consequence of the GlxGr charge 
of ^0 and the electroweak quantum numbers of iJ. As a result, there exists a 
range of parameters in the multi-scalar potential where the vacuum alignment 
of the ^-doublet scalars is essentially independent from the details of the 
VEVs (.^o) and (H). Specifically, we can assume the standard electroweak 
symmetry breaking and allow the field to aquire an arbitrary VEV of the 
order |(Co)| — 10^ GeV. It is therefore sufficient to restrict our considerations 
concerning the vacuum alignment of the ^-doublet scalars to the potential 
V in Eq. (30), in order to determine the range of parameters which leads to 
realistic lepton masses and mixing angles. This analysis will be carried out in 
the following section. 



5 The vacuum alignment mechanism 



We will now determine the vacuum structure emerging after SSB from the po- 
tential V by minimizing each of the individual potentials Vi($, Q) and V2($, Q) 
which appear in Eq. (30). For this purpose it is suitable to parameterize the 
VEVs of the doublets $ and Q as 
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(0a) 1 




^e^"^! cos 
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= 










^^e^'^'i sin a ^ 












^e'<^2 cos 












^e''^2 sin i5 J 


= V2 












S/3 



(31a) 



(31b) 



where ^1,^2 are positive numbers and (/?!, (/p'^, (/92, denote the phases of the 
VEVs. For convenience, we will work with the relative phases Lp = ip[ — (pi and 
ip = ip'2 — ^2- For all potentials Vo($) which appear in Eq. (30) we take the 
quadratic couplings Aq, defined in Eq. (27), to be negative. Furthermore, we 
assume for all potentials V3($, Q) in Eq. (30) the quartic couplings ^3, defined 
in Eq. (29c), to be positive and sufficiently large compared with the remaining 
quartic couplings in Eq. (30). This ensures non- vanishing VEVs f 1, ^2 7^ and 
vacuum stability. Then, the possible physical vacua can be analyzed using the 
parameterization of Eqs. (31) where vi and V2 are kept fixed, i.e., one can 
treat each of the fields $ and as a non-linear cr-model field and restrict the 
analysis to the {a, (p, '0)-parameter-subspace. 



In Eqs. (27) and (29c) we observe that the potentials Vo($) and Vs{^,Q) 
exhibit an accidental ?7(l)^p(,-symmetry where ^/(l)fcc = U{l)a x U{l)p x 
U{l)^p X U{1)^. As a consequence, the potentials Vo{^) and V3{^,n) which 
appear in Eq. (30) have no infiuence on the vacuum ahgnment of the SM 
singlet scalars. Therefore, we can without loss of generality discard them for 
the rest of our discussion. 



Assume first that the fields $ and Q both transform according to the same 
irrep [i = 1,2,3,4). For the most general potential involving only these 
fields we will denote by Va($, ^) the part which breaks the U (l)acc symmetry. 
From Eq. (29a) we find that Va($, ^) can be organized as a sum Va($, ^) — 
Va($, fl) + Vb{^, O,) of two potentials which exphcitly read 



^ More generally, one can view the parameters a, f3, ip, and as the VEVs of scalar 
fields a.{x), P{x),(p{x), and t/jix), respectively. The scalar field a.{x), e.g., is then 
the coordinate of the manifold of cosets U{l)'^^/U{l)a (which is trivial here) at 
each point of spacc-timc. An alignment of the VEVs of $ and O, with respect to a 
happens when in the lowest energy state a{x) provides only a non- Zmear realization 
of the group C^(l)acc (corresponding statements apply to the fields P{x),(p{x), and 
Vi(x)). 
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+d7(0l06 - 0j0a)(t^lt^6 - ^^ftt^a), (32) 

where the coefficients di,...,dj are some real- valued numbers. The poten- 
tial VA($,f2) depends only on the angles a and (3 whereas Vb{^,^) is, in 
addition, also a function of cp and ■0. If the fields satisfy $ e and 
n e {Ql, Ql, Ql, Ql} the analogous C/ (l)j[cp-breaking term yB(^, must have 
dg = dj = 0. In this case, the relative phases (/? and ip are not correlated in 
the lowest energy state. In Eq. (30) we will assume for each of the differ- 
ent t/(l)f[j,p-symmetry breaking parts Va{^,^) that di,d2 < and that the 
condition 

did2 > 4^3, (33a) 

is satisfied. Additionally, we assume that for all possible terms Vb{^,^) in 
Eq. (30) the coefficients d^ and 0^5 are negative and that they also obey the 
constraint 

{-2d4vt + \de\vlvi){-2d5V^ + \de\vlvl) > d'^^vfvi (33b) 

As shown in Appendix B, the conditions formulated in Eqs. (33) enforce the 
nonzero VEVs of the component fields to satisfy the relations 



(U=±(ei6) (^ = 1,2), (34a) 
(gy=±(gj,) (j = l,2), (34b) 

{ql)^±{ql,) (A; = 1,2,3,4), (34c) 

i.e., within each of the scalar ^-doTiblcts the VEVs of the component fields 
are relatively real and exactly degenerate (up to a possible relative sign). 
Note that all terms in the potentials Vb($,J1) which are multiplied by the 
coefficient dr must vanish in the lowest energy state and, hence, cannot con- 
tribute to the minimization of the scalar potential. For the potential V we 
furthermore choose in both of the terms Vb{Qi,Qi) and Vb{^i,^2) fhe corre- 
sponding coefficients to be positive. In contrast to this, the non-vanishing 
coefficients d^ in all of the remaining terms Vb($, Q,) of V are all assumed to 
be negative. Then, the absolute minimum of the multi scalar potential is in 
addition to Eqs. (34) furthermore characterized by the relations 



We consider here only the tree-level approximation. 
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(92a) _ Mil = Mil = Mil 

(qI,) {qI,} {qD {qD' 



(35a) 
(35b) 
(35c) 



i.e., the relative orientation of the VEVs of the component fields within a 
specific doublet is equal for all doublets transforming under the irrep D4 and 
opposite for the pairs of doublets transforming under or D^. 
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We suppose that the VEVs of the fields ^0, Ci, and {2, which are responsible for 
the generation of the neutrino masses, are all of the order of the electroweak 
scale 

|(eo)|^|(UI=^l(6a)|=^10^GcV. (36) 

In contrast to this, we assume that the link fields of the deconstructed extra- 
dimensional gauge symmetries and all aquire VEVs of the same order 
of magnitude at some high mass scale somewhat below Mf and thereby give 
rise to a small expansion parameter 



A ~ ~ ~ 0.22, (37) 

where j = 1,2, and k — 1,2,3,4. Small hierarchies of this type can emerge 
from large hierarchies in supersymmetric theories when the scalar fields aquire 
their VEVs along a D-flat direction [58]. Note in Eq. (37) that A is given by 
the Wolfenstein parameter [59] which approximately describes the mass ratios 
and CKM mixing angles in the down-quark sector [60] as well as the mass 
ratios in the charged lepton sector [61].^^ The mass and mixing parameters 
of the charged leptons are determined in the next section. 



6 The charged lepton rticiss matrix 



Consider the Yukawa interactions of the charged leptons 

= L:H^i^Ef3 + h.c., (38) 



The terms associated with the coefficients de actually represent a spin-spin- 
interaction in a version of the Ising-model, known from ferromagnetism. The topol- 
ogy here is of course unfamiHar since all "spins" couple with equal strength. 

An Ansatz where the Wolfenstein parameter is also used to describe neutrino 
mixing and leptogenesis has recently been presented in Ref. [62] . 
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E„ 


E^ 




Gt X Gn 










(-2,2) 


(-4,4) 


(-2,3) 


(-2,3) 


l; 


(-1,0) 


(-3,2) 


(-1,1) 


(-1,1) 




(-1,0) 


(-3,2) 


(-1,1) 


(-1,1) 



Table 3 

Gl X Gr charge structure of the charged lepton-antilepton pairs. 

where denotes an effective Yukawa operator and a, /3 = e, ji, t. The Gl x 
Gr charge structure of the charged lepton-antilepton pairs is shown in Table 
3. One of the possible lowest-dimensional contributions to the operator is 
given by 



(39) 



since it yields an invariant under both Gl x Gr and ^. Furthermore, invari- 
ance under application of the operator T>[a4) requires ff^^ = i^^^ = 0. Since 
the transformation T>{a2) permutes <-> ^^^t ^^'^ ^fie ^ ^^re ^ig^ 
flip is possible due to .03(02) = D{Gb)) it follows that = = 0, i.e., in 
the flrst row and column of the charged lepton mass matrix only the (1,1)- 
element is non- vanishing. From ^7(1) gauge- invariance we have seen in Sec. 3.1 
that the lowest-dimensional mass terms in the //-r-subsector of the charged 
leptons are given by the dimension-six and dimension-eight operators shown 
in Fig. 4. In theory space, the effective Yukawa couplings of the dimension-six 
and dimension-eight operators are identified with the Wilson loops around the 
plaquettes associated with the deconstructed extra-dimensional gauge symme- 
tries G^ and (see Figs. 5 and 6). The generahzed Wigner-Eckart theorem 
imphes that each of these Wilson loops corresponds in ^-space to a set of 
irreducible Yukawa tensor operators spanned by the irreps 2^ and 2^. These 
operators can be quickly determined by first noting that under ©(04) the ef- 
fective couplings ff^^ and undergo a sign flip whereas all link flelds of 
n^^j^Gj^ and II^^j^G^ transform trivially. As a result, we have = — 0, 
i.e., the charged lepton mass matrix is diagonal. Then, testing for invariance 
under the action of 01,02,03 G ^ gives that each set of irreducible Yukawa 
tensor operators transforms according to a representation D5 of ^ which is in 
matrix-form defined by the generators 




DM = 




0)5(03) = 05(04) = I2. (40) 



To leading order, the generators in Eq. (40) act on five independent doublets 
(V'a i^h^y product functions ip^^^ and where i = 0, 1, . . . , 4, which form 
the basis of five distinct carrier spaces \/(*). Here, the $f-doublets (^(^\ ^J^V 
correspond for i = and i > 1 to the Wilson loops around the plaquettes 
associated with G^ and G^, respectively. The basis functions for the different 
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^b 


T/(0) 
V 




Qlb<l2b 


T/(l) 




llblibllblib 


F(2) 


llaliallb^lb 


1lbl2bl3alla 


y(3) 


(iallbliaQlb 


1lbl2al3bl4a 




^Ial2bl3blla 


1lbl2a'l3al4b 



Table 4 

The basis functions V'a'* and ijj^^ {i = 0, 1, ... ,4) of the leading-order set of irre- 
ducible Yukawa tensor operators (^|;), where k,l = /i, r, in terms of the component 
fields of Qj = {ql, qlY {i = 1, 2) and Qf = q^,f {i = 1, 2, 3, 4). 



carrier spaces of are given in table 4. The allowed types of basis functions in 
table 4 follow from invariance under application of the transformations X>(a3) 
and V{ai). First, application of 2^(03) yields that in each basis function the 
number of the a-components is even. Second, for given i = 0, 1, . . . , 4 invari- 
ance under the transformation V{ai) requires the index structure of the basis 



functions ijj^'^ and t/j^ to be of such a form that tjj^'^ and i/j^ get interchanged 
when the indices a and b are permuted. Then, expressed in terms of the sets 
of irreducible Yukawa tensor operators, the effective Yukawa coupling matrix 
in the //-r-subsector of the charged leptons reads 







1 



-1 



+1 I ^, 
1/ I 1 



(41a) 



where k,l = ij,,t and Y^g for i = 0,1,..., 4 denotes the effective Yukawa 
couplings 

Vw/(M/)2 (i = 0) 



eflf 



y«/(My)4 (i = l,2,3,4) 



(41b) 



where Y^^^ (i = 0, 1, . . . , 4) are order unity Yukawa couplings. The 2x2 diago- 
nal matrices in Eq. (41a) summarizing symmetry-related geometric factors, are 
the Clebsch-Gordan coefficients of the effective Yukawa coupling matrix {^ii). 
Furthermore, the effective Yukawa couplings Y^^ ^ characterized by the outer 
multiplicity label i, are the reduced matrix elements of the Clebsch-Gordan 
coefficients and parameterize further information about the physics at the 
fundamental scale Mf. Prom Eqs. (34) we find that after SSB the vacuum 
alignment mechanism of Sec. 5 ensures that the VEVs of the basis functions 
in table 4 are - up to a possible relative sign - pairwise exactly degenerate 



(42a) 
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where i = 0, 1, . . . , 4. In addition, Eqs. (35) relate the orientations of the VEVs 

by 



(0)\ 



(42b) 



where i = 1,2, 3, 4. Substituting Eqs. (42) into Eqs. (41) we observe that after 
SSB the set of irreducible Yukawa tensor operators {^li) can take one of the 
following two forms 



y(0);,2 




i=l 



1 

,0 1, 




(43) 



where k,l — ijl,t and the expansion parameter A ~ 0.22 of Eq. (37) has been 
used. In Eq. (43) it is important to note that the Clebsch-Gordan coefficients 
are added or subtracted depending on their outer multiplicity label: if the 
Clebsch-Gordan coefficients are added (subtracted) for i = then they are 
necessarily subtracted (added) for i = 1,2,3,4. As a result, Eq. (43) shows 
that the vacuum alignment mechanism generates a hierarchical pattern in the 
/i-r-subsector of the charged leptons via a cancellation of some of the Clebsch- 
Gordan coefficients in the lowest energy state. For definiteness, let us choose 
in Eq. (43) the solution with the signs "+" for i = and "— " for i = 1,2, 3, 4. 
Taking everything into account, the full leading order charged lepton mass 
matrix M.^ emerging after SSB is given by 



Mi 



0^ 
A2 
1 



(44) 



where rrir is the tau mass and only the orders of magnitude of the matrix 
elements have been indicated. The masses and the mixing of the charged 
leptons can be calculated by diagonalizing Aie-M-l- Denoting the electron and 
muon masses by rUe and m^, respectively, the mass spectrum described by 
is found to be 

^ ~ A^ ^ ~ A2, (45) 
rrir rrir 

which approximately fits the experimentally observed values [61]. As for the 
mixing angles of the charged leptons practically vanish, the experimentally ob- 
served large leptonic mixing must stem from the neutrino sector. The neutrino 
mass and mixing parameters will be determined in the next section. 
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Fig. 7. Moose diagram for the location of scalars and chiral fermions in the orbifold 
extra dimensions associated with the gauge groups and G^. 

7 The neutrino mass matrix 



7.1 Aliphatic models 

So far, we have examined the connection between dynamically generated extra 
dimensions compactified on and the hierarchical Yukawa coupling matrix of 
the charged leptons. In such a geometric approach to small Yukawa couplings 
it is generally interesting to relate the leptonic mass and mixing parameters 
to different topologies in theory space. For this purpose we will now, following 
Sec. 2.2, consider for each of the gauge groups = n^\Gf and G'^ = U^JiGf 
the aliphatic model for fermions in the latticized S^/Z2 orbifold extra dimen- 
sions. Since these fermions are SM singlets, we can identify them with "right- 
handed" {i.e., SU{2)l singlet) neutrinos. We suppose that the fundamental 
Froggatt-Nielsen states of the charged lepton sector transform only trivially 
under G^ and (see Sec. 3.1). Then, the orbifold extra dimensions are (at 
tree-level) completely decoupled from the charged leptons and can only be 
experienced by the neutrinos. 

At this stage, we allow the number of sites A^3 of the gauge group G^ to be 
large but leave it yet unspecified^^. For the gauge group G^, however, we 
assume a very coarse latticization where the aliphatic "chain" of consists 
only of the N4 — 2 sites Gf and Gf. By assigning ^0 and ^1 the Gf charge +1 
these fields are put on the 1st site associated with G^. Additionally, we assign 
,^2 the Gf charge +1 which locates the field at the site representing Gf. The 
corresponding moose diagrams are shown in Fig. 7. Since the charged leptons 
cannot experience these topologies, the site-fields io,ii, and ^2 are now only 
relevant for the generation of the neutrino mass matrix and we can discard 
them in the discussion of the charged lepton mass matrix. 



In Sec. 7.4 we will show that the number of sites parameterizes the solar 
neutrino mass squared difference Amg. 
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1.2 The one- generation- case 



Let us first restrict to the case of one generation by considering the coupUng 
of the ^f-singlet neutrino to the fermionic site variables of the orbifold extra 
dimensions. From tabic 1 wc conclude that x gauge-invariancc requires 
the relevant effective Yukawa interaction of the electron lepton doublet 1^ with 
the bulk-fermions to be of the type 



(46) 



where Y\ and are order unity coefficients, H is the charge-conjugated Higgs- 
field H = ia^H*, and Mip 10^^ GeV denotes the seesaw scale. Suppressing 
the index j = 3, the Lagrangian for the bulk and brane fields is given by 
^1 + ^2, where =^1 has been defined in Eq. (10). As a result, the A^3 x A^3 
mass matrix which emerges from jSfi + ^2 after SSB is given by 



1 



IfL *1L *2L *3L • 

hi 

MiF /i2 

Mf -Mf 
M/ -Mf 



\ 



(47) 



where we have introduced the VEVs hi = Yi{H*), /12 = ^2(^0); and use 
/ii ~ — 10^ GeV, i.e., the associated mass terms are generated at the 
electroweak scale. Note that the mass matrix in Eq. (47) has been displayed 
in a basis where the active neutrinos have right-handed chirality. Integrating 
out the heavy ^-singlet Ip, yields the effective {N^ — 1) x {N^ — 1) mass matrix 



••• 

Mf -Mf • • • 
Mf -Mf ■ ■ ■ 



V 



(48) 



7 



where rrij, = hih2/Mip ~ 10^^ eV denotes the absolute neutrino mass scale. 
The mass m,y lifts the zero eigenvalue in the KK mass spectrum of the left- 
handed fields to a small but non- vanishing value which can be determined by 
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diagonalizing the (A^s — 1) x (A^3 — 1) matrix 



ml rriuMf 
m^Mf 2Mf -Mf ■ ■ ■ 
-M? 2M? ••• 



M2MI 



(49) 



v 



■7 



From Eq. (49) it is readily seen that the mixing of Ij with the right-handed 
fields "i/nR is described by angles ~ m^/Mf which vanish in the limit Mf — > 00. 
For <C Mf the non-zero heavy masses are approximately given by the 
KK spectrum in Eq. (15) where n — 1, . . . , — 2. The lightest eigenvalue 
(mo)^ of can be determined by integrating out the invertible heavy 

(iVs — 2) X [N^ — 2) submatrix in the down-right corner ofM2Ml in Eq. (49). 
Taking into account that the (1, l)-element of the inverse of this submatrix 
is equal to ^^M^^, we obtain in the limit rrii, <C Mf for the lightest mass 
eigenvalue 



where the second equation matches onto the continuum limit. In the classi- 



cal 5D continuum theory the mass (or volume) suppression factor ~ ^jMfR 

emerges from the normalization of the wave-function of the right-handed neu- 
trino propagating in the bulk [63]. In our deconstruction setup, Eq. (50) states 
that for a given number of sites there is a non-decoupling of the deconstruction 
physics from the low-energy theory. Since, as already stated after Eq. (49), the 
mixing of the active SM lepton 1| with the right-handed fields \E'„^ vanishes 
in the limit Mf — 00, we can repeat the above construction for additional 
aliphatic models which, after integrating out the corresponding KK towers, 
give rise to further light Dirac neutrino masses. In this way, one can build up 
a Dirac neutrino mass matrix, which reflects the properties of the dynamically 
generated extra dimensions. 

7.3 Inclusion of the 2nd and 3rd generations 

The full neutrino mass matrix emerges by inclusion of the SU{2) lepton dou- 
blets of the 2nd and 3rd generation, which are combined into the ^-doublet 
2(. Actually, gauge- invariance under Gl x Gr allows for the flelds ^2 and ^3 
Yukawa interactions with 2^ which read 



^3 = Y^2'iH*2FL + 142^^6 + ^52fr6*il + M2f2fr2fl + h.c, (51) 



where Y3, 14, and Y5 denote order unity Yukawa couplings, (^1) ~ (^2) — 
10^ GeV, and M2F — 10^^ GeV. In the construction of Sec. 2.2 we associate 




(50) 
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with the Weyl spinor a two-site model for where resides on the 
site corresponding to Gf (sec Fig. 7). In complete analogy to the calculation of 
the light mass mo in Sec. 7.2, one can in the combined system {J2i=3,A-^i) + 
^2 + -^3: where J2i=3,4^i summarizes the ahphatic chains of G^ and G^, 
integrate out the heavy vectorlike degrees of freedom. As a consequence, in 
the basis where the VEVs of and ^2 are described by Eqs. (34a) and (35a) 
the resulting 3x2 light Dirac neutrino mass matrix can be written as 



Mn = 



e 1 



(52) 



where the neutrino expansion parameter e ~ 1/ \/iV3 — 1 and the order unity 
Yukawa coupling p = 0{1) are both real quantities and m,^ denotes the ab- 
solute neutrino mass scale. In Eq. (52) all phases have been absorbed into 
the right-handed neutrino and charged lepton sectors implying the practi- 
cal absence of CV violation in neutrino oscillations. The Dirac neutrino mass 
matrix in Eq. (52) has the important property that within each column 
the flavor symmetry $f enforces the 2nd and 3rd elements to be relatively real 
and exactly degenerate in their magnitudes. Thus, M.d describes an exactly 
maximal i/^-i/^-mixing. The structure of Ain is familiar from models of sin- 
gle right-handed neutrino dominance [28] which provide an understanding of 
normal hierarchical neutrino mass spectra in the context of the MSW LMA 
solution. 



7.4 Neutrino masses and mixing angles 



The neutrino masses and leptonic mixing angles are determined from Eq. (52) 
by diagonalizing the matrix 



-pe- 



\ 



pe^ 1 + 6^1- £2 

-pe^ 1 - eM + 



(53) 



The matrix in Eq. (53) is brought on diagonal form by a rotation of the active 
neutrino fields in the 2-3-plane through an angle 6*23 = 7r/4 followed by a 



This is possible since the charged lepton mass matrix Ai^ in Eq. (44) is of diagonal 
form. 



28 



rotation in the 1-2-plane through an angle 



9i2 = arctan 



(2x/2)-^ (p2_2 + ^(2-p)2 + 8) 



(54) 



Hence, the reactor mixing angle ^13 exactly vanishes in agreement with the 
CHOOZ reactor neutrino data which sets the upper bound j^ia] < 9.2° [65]. 
The neutrino masses exhibit the normal hierarchy 



mi = 0, 1712 = mve\j2 + p^, 777,3 = V^rn^^ (55) 
which gives for the solar and atmospheric neutrino mass squared differences 

A777I = 777 ^62(2 + p2), A7772^^ = 2777^ - A777|. (56) 

Using the upper bound A777q < 1.9 x 10~^ eV^ [17] and the best-fit value 
A7772,^ ~ 2.5 X 10--'^ cV^ [9] we obtain 777j^ ~ 0.04 eV which is consistent 
with the value Mp^ — — 10^^ GeV for the seesaw scale. Without tuning 
of parameters we have p — 1 and e = 1/ 1/^3 — 1 which gives for the solar 
neutrino parameters the values 

3 Att?^ 1 
^<^2N;^' ^12 = arctan ^ 35°, (57) 

where we have used in the first equation the hierarchy A777q ^ A777^jj^. At 
3(7, the combined solar and KamLAND neutrino data allows for A777q the two 
regions 1^ 5.1xlO-^eV2 < A777I < 9.7xlO-^eV^ (LMA-1) and 1.2x10"^^ < 
A777I < 1.9 X lO-'^eV^ (LMA-II) [17]. Matching onto these values requires 

A?^3 = 57±17 LMA-I, Ar3 = 27±6 LMA-II, (58) 

where we have set in Eq. (57) the atmospheric mass squared difference equal 

to the best-fit value A777^^j^ = 2.5 x 10~^ eV^. In short, the presently allowed 
ratios A777q/A777^^.jjj implied by the LMA-I and LMA-II solutions already sig- 
nificantly discriminate between the associated radii N^/ {271 gv) of the dynami- 
cally generated /Z2 orbifold. At this level, the neutrino expansion parameter 
is 0.12 < e < 0.16 (LMA-I) or 0.18 < e < 0.23 (LMA-II), which is comparable 
with the Wolfenstein parameter A ~ 0.22 of the charged fermion sector. For 
the non-fine-tuned solar mixing angle 612 = arctan l/-\/2 in Eq. (57) we find 
from the analysis in Ref. [18] that a number of 

55 ± 8 LMA-I (@ 90% C.L.) (59) 

lattice sites yields the MSW LMA-I solution within the 90% confidence level 
region. In general, for both the LMA-I and the LMA-II solution, the dynamical 



At tree- level, sub-leading corrections may come from the charged lepton sector. 
We adopt here the nomenclature of Ref. [16]. 
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generation of the solar mass squared difference Am? via deconstruction in a 
flat background requires a relatively fine-grained latticization of the associated 
/Z2 orbifold with roughly 10^ — 10^ lattice sites. 



8 Summary and Conclusions 



In conclusion, we have presented a model for lepton masses and mixing an- 
gles based on a non-Abelian discrete flavor symmetry and U{1) charges in a 
deconstructed setup. The model applies the vacuum alignment mechanism of 
previous models [33, 34] in order to predict an exactly maximal z/^-z/T--mixing 
as well as the strict hierarchy -C between the muon and the tau mass. 
The non-Abelian discrete symmetry is identified as a split extension of the 
Klein group Z2 x Z2 where the lift of every fermion and scalar representation 
is equivalent with the dihedral group ^4 of order eight. The charged lepton 
masses are generated by the Froggatt-Nielsen mechanism which is given a 
geometric interpretation in terms of deconstructed or latticized extra dimen- 
sions compactified on the circle . In theory space, the muon and tau masses 
correspond to Wilson loops around the plaquettes associated with the decon- 
structed extra-dimensional U{1) gauge symmetries. As a result, the model 
gives the realistic charged lepton mass ratios irie/mT- ~ and m^j_/mr — A^, 
where A ~ 0.22 is the Wolfenstein parameter. Since the charged lepton mass 
matrix is of diagonal form, the leptonic mixing angles stem entirely from the 
neutrino sector. Enforced by the symmetries, the vacuum structure yields an 
exactly maximal atmospheric mixing angle Qiz — t^I^ and a vanishing reactor 
angle ^13 = 0. In addition, all CP violation phases vanish due to the symme- 
tries. The model provides an order of magnitude understanding of the solar 
mixing angle 6*12, which is predicted to be large but not necessarily close to 
maximal. Specifically, without tuning of parameters (i.e., by choosing univer- 
sal values for all real Yukawa couplings of order unity) the model yields the 
solar mixing angle Q\2 = arctan l/\/2 ~ 35°. The neutrinos exhibit a normal 
mass hierarchy through single right-handed neutrino dominance which is real- 
ized by the propagation of right-handed neutrinos in latticized S^/Z2 orbifold 
extra dimensions. Here, the solar mass squared difference Am^ is suppressed 
against the atmospheric mass squared difference Am^^^ by the discretized 
analogue of the volume factor, known from the classical theory. For a latti- 
cization of the orbifold with 10^ — 10^ lattice sites the model yields without 
tuning of parameters the MSW LMA solution (LMA-1) of the solar neutrino 
problem within the 90% confidence level region. 
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Fig. A.l. Horizontal plane with the system of two-fold axes for ^4. 
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A The dihedral group ^4 

The dihedral groups where n = 2, 3, . . ., are the point-symmetry groups 
with an n-fold axis^^ and a system of 2-fold axes at right angle to it. The 
group is therefore the symmetry group of a regular n-gon. These groups 
contain 2n elements and for n > 2 they are non-Abehan (^2 is isomorphic 
with the Klein group x Z^). In Fig. A.l the horizontal plane is shown for 
the case n = 4, where a, a', 6, and h' denote the four two-fold axes and the 
4-fold axis is perpendicular to the paper. We denote by C„ the operation of 
rotation through 27r/n about the principal axis. The fc-fold application of this 
transformation will be written as C\ and the identity transformation C" as 
E. The rotations through tt about the axes a, a', 6, and h' will be referred to 
as Ca,Ca',Cb, and Cb', respectively. Then, the dihedral group ^4 has eight 
elements in following five classes: 

E; C4,Cl; Cj; Ca,Cb; Ca',Cb'. (A.l) 

Adopting the notation of Ref. [66] we will refer to the five classes which are 
associated with the sequence in Eq. (A.l) as 6*4(2), C|, C2(2), and C2'(2). 



This is also called the principal axis. 
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Table A.l 

Character table for the group ^4. 
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Table A.2 

Multiplication table for the group ^4. 



Calling the 4 irreducible singlet representations 1a, 1b, Ic, and Id respec- 
tively, the decomposition of the product of the two-dimensional irrep 2 reads 



2 X 2 = 1a + 1b + Ic + Id- (A.2) 

The character table for the group ^4 is given in table A.l. Denoting 2 as (a, h) 
we have for the singlet representations 



lA = aia2 + feife2, (A.3a) 

1b = 0162 - , (A. 3b) 

Ic = aia2 - &ife2, (A. 3c) 

Id = 0162 + 02^1 • (A. 3d) 



Prom the character table of ^4 one determines the decomposition of the prod- 
uct of any two representations as shown in table A.2. In the vector represen- 
tation of ^4, the 2x2 representation matrices corresponding to the different 
classes can be written as 
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E 



C4(2) 



C2(2) 

C2'(2) 




D{Cl) 



D{C\ 



-1 

1 




(A.4) 



Before concluding this section, let us construct the dihedral groups from semi- 
direct products. For this purpose, let A/" ~ Z„ for any n e N, let ~ Z2, 
and let the map (p : H ^ Aut(A^) send /i e if to the automorphism in 
Aut(iV) sending each element of A^ to its inverse. Then, the external semi- 
direct product N xi^ H oi N hj H with respect to (p is the dihedral group 
One can also define the infinite dihedral group = N H, where N Z 
and the group H and the map are as above. 



B Minimization of the tree-level potential 



We will rewrite the potential Va($, ^) in Eqs. (32) in terms of the parameter- 
ization in Eqs. (31) as follows 



Va = d.vtclsl + d^vtclsl + d^vlvlicl - sl){cl - s^), (B.l) 
where Sa = sin(Q;) and = cos(q;) (correspondingly for /3). Hence, it is 



— ^2divf{c'^^Sa - CaS'^a) " 4:d3vlv2CaSa{c^p - S^), 
2d2V2{CpS/3 - CjiS^) - 4:d3vlv2Cf3S(3{cl - sl). 



dp 



As a result, at (a,/5) = (f , f ) it is {da,dfi)VA = (0,0), i.e., (a,/3) = (f , f ) is 
an extremum of the potential Va- For the second derivatives it follows 
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= 24^^(4 + 4 - 644) - ^d^^l^iic} - 4) (4 - 4), 



d'VA 

da dp 



16d3vlv2CaSaC(3Si3. 



At {a, (3) = (f , f) we therefore obtain for the matrix of second derivatives of 
the potential 



da 8(5 



-2divf 4(^344 



\^4(i3?;^f| —2d2V2 



(B.2) 



where, due to the choice di,d2 < 0, the diagonal elements are positive. Hence, 
if the parameters obey the condition 

did2 > 4:dl, (B.3) 

the matrix in Eq. (B.2) is positive definite, i.e., the modes which oscillate in 
the (a, /3)-subspace around {a, P) — (|, |) have positive masses and {a, P) — 
(j, j) is a minimTim of the scalar potential V^. Let us now rewrite the part of 
the multi-scalar potential Vb{^, fl) in Eqs. (32) using the parameterization in 
Eqs. (31) as follows 



Vb — 2(i4f ^4'5aCos(2(^) + 2(^544=^3 cos{2tjj) 

+4:dQvlv2CaSaCf3S 13 COs{'^) COs(l/)) 

-Ad'jvlvlcaSaCpSp sm^if) sin(z/;), (B.4) 
where we have used the notation of Eq. (B.l). Hence, one concludes 



= 4(i4t^i(4sa - Cc,s%) cos(2v9) + 4:d(ivlvl{cl - s^)c/3S/3 cos{ip) cos(t/') 

-4.divlvli.cl - 4)c/3S/3 sin((^) sin(V^), 
= 444(4^/3 ~ C0^V) cos(2V') + Ad(yvlvl{c^i3 - s'p)CaSa cos((/?) cos(V') 



da 

-Mivfviicl - sDcpSfj sin((^) sin(V^), 

dp 

-MivfvKcl - sl)cc,Sa sin{ip) sin(t/'). 



and 
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= -4d4V^clslsm{2(p) - Ad^vlvlcaSaCpSp sin(v3) cos(V') 

-Ad-jvlvlcaSaCpSp cos((^) sin(^), 
= -Adr,v^c^pSpSm.{2'ip) - AdevlvlcaSaCisSp cos{(p) sm{ip) 

—Adivlvlca-SaCpSp sm.{ip) cos(-(/'). 

As a result, at the points (ct, (3) — (|, |), where <f,i>& {0, tt}, it is 

(a«,a;3,a^,a^)ys = (0,0,0,0), 

i.e., these points are extrema of Vb- Furthermore, one finds at (a, /3) = (f , f ) 
vanishing mixed second derivatives 

d^VB d'VB d'VB O^Vb ^ Q 



d(pda dip da dLpd(3 dip 8(3 

implying that the matrix of the second derivatives of Vb with respect to the 
parameters a, /9, ip breaks up into a block-diagonal form with submatrices 
which respectively correspond to the subspaces (a,/?) and {(p,ip). The second 
derivatives of Vb with respect to a and /3 are 



d'VB 



d-'VB 



d'VB 

dadp 



Ad4V^{sl + 4 - 6clsl) cos(2(/7) - l^d^vlvlcaSaCpSp cos((/7) cos(V') 
-\-lQd^vlv2CaSaC0S|3 sin((^) sin(-?/'), 

4^5^2(5^ + dp — 6c|s^) COS(2'0) — 16devfv2CaSaCf3Sp cos((^) cos('0) 
+ 16dYvlv2CaSaCi3S 13 sm{(p) sin('0), 

4rf6?;it;2(c« - sl){cl - si) cos{ip) cos(^) 
-Adrvlvlicl - sl){cl - s^) sm{ip) sm{ip). 

(f'f)' where (p,ip E {0,7r}, the matrix of 



Therefore, at the points {a, (3) 
the second order derivatives is 



' O^Vb ' 
da 8(3 



—d^V2 — adQvfv2^ 



(B.5) 



where a = cos((/9) cos{ip) = ±1 can take either sign for ip^ip E {0, tt}. However, 
from Eq. (B.4) it is seen that the product d^a must be negative in the lowest 
energy state, i.e., the sign of determines whether 93 = ip+k-2T{ oitp = ip+k-n 
for some integer k. The matrix of second order derivatives can therefore be 
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rewritten as 



9Vi 



B 



da d(3 



\ 



-d4vf + \dQ\vjv2 ^ 

—d^V2 + \dQ\v\v2j 



(B.6) 



where ^4,(^5 < 0, i.e., the matrix is positive definite. The second derivatives 
of Vb with respect to </? and are 



B 



B 



di/jdcp 



-Sd^v^cl^sl^ cos(2(^) 

—AdQvfvlcaSaCpSp cos{(f) cos{ip) + AdrvlvlcaSaCpsp sm{(f) sin('0), 

— 8(i5f 2C^s^ cos{2ijj) 

—AdQv\v\caS aCjjS 13 cos((^) cos('0) + Adiv^vl^CaSaCps jj sin((^) sin('0), 

A:dQv\v2CaSaCl3Si3 Sin((/7) sin('0) — A:d'jv\v2CaSaCl3Sp C0s(</7) COs('0). 



At the points = (f , f ), where (p,ip E {0,7r}, the matrix of the second 

derivatives of Vb with respect to ip and ip reads 



\dipdijj J 



—2d4vf + \dQ\vlv2 ^dful'ul 

±d7vfv2 — 2o?5f| + \a\devlv2 



(B.7) 



where ^4, ^5 < 0, i.e., the diagonal elements are positive. In Eq. (B.7) we have 
already used that the potential is minimized when d^cr is negative. 

Taking everything into account, if the coefficients in the multi-scalar potential 
satisfy besides Eq. (B.3) also the condition 



—2d4Vi + \dQ\v^V2){—2d5V2 + \dQ\v^V2) > d.^ViV2, 



(B.8) 



then the matrix in Eq. (B.7) is positive definite, i.e., all modes oscillating in the 

(a, /9, (/?, ■?/')-subspace around the points (a,/9) = (f-f), where (/?,-?/' G {0,7r}, 
have positive masses and hence these points are indeed local minima of both 
the potentials Va and Vb, i.e., they locally minimize the term Va($, ^) which 
breaks the accidental C/(l)^[^.j,-symmetry. 
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